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ABSTRACT 


An experimental investigation was conducted to examine 
E Convective heat transfer ln straight and curved ducts of 
rectangular cross-section. The experimental configuration 
Was modeled as infinite parallel plates with one wall at a 
constant temperature and the opposite wall adiabatic. 

The experiments were conducted at steady state for tur- 
bulent flow. Average Nusselt numbers were used to compare 
the heat transfer characteristics of the straight and curved 
Sections. The development of Taylor-Gortler vortices in the 
Curved section was shown to enhance the heat transfer rate in 
the curved section as compared to that of the straight section 
approximately 15 to 20 percent. Improved heat exchanger 
EEEE ns and improved cooling of turbine blades are two 415- 
ciplines that could benefit from a better understanding of 


the effects of curvature on the rate of heat transfer. 
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I ZENTERODUCETTON 


EN AYLOR-GORTLER VORTICES 

since the early part of the twentieth century, consider- 
able research has shown that fully developed laminar flow 
along a concave wall does not remain two-dimensional [Refs. 

1, 2, 5]. The flow instead forms a system of spiral vortices, 
Of counter rotating pairs, whose axes are in the principle 
flow direction. This phenomenon is the result of the varia- 
ems inthe centrifugal forces acting on the fluid particles, 
ARG is Known as Taylor-Gortler vortices. Figures 1 and 2 
MS trate the type of fluid motion just described. 

In a channel that is curved in the streamwise direction, 
eee fluid particles located in the center of the flow cross- 
Section are subjected to higher centrifugal forces than those 
iud particles traveling along the channel's boundary wall. 
As a result, the fluid in the center of the channel moves out- 
wardly towards the concave boundary. As the process continues, 
the fluid particles near the boundary wall, move in a spanwise 
ection, and finally radially inward replacing the outwardly 
moving particles. These particles then come under the same 
нив Роса! forces and the process repeats itself. The re- 
sulting cyclic motion, causes the formation, and propagation, 
Berne counter rotating Taylor-Gortler vortices which are con- 


sidered primarily a laminar flow phenomenon. 
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It has been observed [Ref. 4] that the heat transfer rate 
from flow along a concave curved wall is greater than that 
TS flow along a straight wall. It is thought that the addi- 
тај Mixing provided by the secondary motion of the Taylor- 
Eder vortices may account for this. It is also thought that 
the cross-hatching that has been observed in reentry vehicles 
can be explained in part by the presence of streamwise vortices 
similar to the Taylor-Gortler vortice phenomenon. 

Ге аге many possible applications that could result from 
a more thorough understanding of Taylor-Gortler vortices and 
their effect on heat transfer and fluid flow characteristics. 
By taking advantage of the increased heat transfer rate, im- 
proved heat exchanger designs and improved turbine blade cool- 


EuSEcould result [Refs. 5, 6, 7]. 


ПИ HISTORY 

The instability of an inviscid fluid flowing past a curved 
boundary, was first considered by Lord Rayleigh in 1916 [Ref. 
8]. By assuming that the fluid was non-viscous, he determined 
ШЕ Тот the motion of an inviscid fluid to remain stable, its 
Emulation must increase with increasing radius. G. I. Taylor, 
925, [Refs. 1, 9], continued this study with an extensive 
analytical and experimental study of viscous fluids. His in- 
Eu Eations focused on the flow between two cylinders. For 
the case where the inner cylinder was rotated while the outer 


cylinder was held stationary, Taylor ascertained that such 
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flows become unstable when the value of a dimensionless para- 
Meter exceeded a critical value of 41.3. The parameter, 


known as the Taylor number is defined as: 


Ta = Re P. 
l 


where 'd' is the width of the gap, assumed small when compared 
EN DU the radius of the inner cylinder, and 'Re' is the 
Reynolds number. The Reynolds number is based on the peri- 
pheral velocity of the inner cylinder. Taylor determined 

ЕШ Гот those cases in which the value of the Taylor number 
exceeded 41.5, a secondary motion developed and the Taylor 
vortices formed. 

Instability of a similar nature is also observed when a 
viscous fluid flows in a curved channel due to a pressure 
gradient acting along the channel wall. This problem was first 
considered analytically by W. R. Dean [Ref. 10] in 1928, for a 
Channel formed by two concentric cylinders, where the radius of 
the innder cylinder was large in comparison to the spacing be- 
tween the inner and outer cylinder walls. Dean concluded that 
there would be an initiation of the flow instability, and the 
puopacation of vortices, when а similar dimensionless parameter, 
the Dean number, exceeded a value of 36. The Dean number is 


defined as: 


2; 





memes d' represents the channel half-width, 'Ri' is the inner 
cylinder radius, and 'Re' is the Reynolds number based on the 
mean velocity of the undisturbed flow. The analytical work 

of Dean was later verified by W. H. Reid [Ref. 11], using an 
approximate solution. 

In 1940, H. Gortler studied the stability of laminar bound- 
ary layer profiles on curved walls under the influence of small 
disturbances. He found that these disturbances were similar 
meme vortices studied by G. I. Taylor. By approximate numer- 
NUEcalculations, Gortler concluded that the disturbances or 
vortices were produced only on the concave boundary walls. He 
So concluded that the overall flow profile appeared to remain 
laminar in nature [Ref. 2]. These results were verified, with 
an exact solution, by G. Hammerlin, and reported by H. Schlichting 
IRef. 12] in 1955. А. М. О. Smith completed an even more ex- 
tensive numerical analysis that further substantiated these 
findings [Ref. 3]. Recently with the use of hot wire anemometry, 
laser doppler systems, and flow visualization techniques, [Ref. 
13, 14], the result of these numerical solutions have been 
demonstrated experimentally. In 1976, Y. Aihara [Ref. 15] con- 
ducted a non-linear analysis of Gortler vortices. 

AS interest began to develop concerning the effects of the 
secondary flows associated with the Taylor-Gortler vortices, 
many studies were published concerning the influences of these 
vortices on the transfer of heat in the laminar and turbulent 


flow regimes. In 1955, F. Kreith [Ref. 16], studied the 
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Eu nce of heat transfer with respect to the curvature of 
the boundary wall for fully turbulent flows. He concluded 
that the heat transfer from the heated concave boundary wall 
Was considerably higher than that transferred from the convex 
boundary wall of the same curvature and under similar turbu- 
Bent conditions. 

In 1965, L. Persen [Ref. 17], considered the special cases 
of very high and very low Prandtl number fluids. He related 
B EMDErease in heat transfer from a curved wall to the pre- 
sence of the Taylor-Gortler vortices. 

There has been only a limited amount of published litera- 
ire dealing with the flow and heat transfer in curved channels 
Of rectangular cross-section and large aspect ratios. Much 
of what has been published involves the development of numer- 
ШЕЕ approximations and solutions for heat and mass transfer 
in curved ducts of various geometries. К. Cheng and M. Akiyama 
[Ref. 18] developed a numerical solution for forced convection 
Meat transfer with laminar flows in curved channels of rec- 
tangular cross-section. However, they were concerned with 
cross-sections with small aspect ratios. 

о, А. А. Shibani and M. N. Ozisik [Ref. 19] solved 
the heat transfer problem between parallel plates with tur- 
bulent flow, for the case of uniform wall temperature. They 
used matched asymptotic expansion techniques for a wide range 


of Prandtl numbers. 
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Other experimental and analytical studies that are worth 
Serre with regard to this present study follow. Y. Mori [Ref. 
20], obtained results for hydrodynamically fully developed 
flows with constant wall heat flux, in curved channels of 
square cross-section. WW. M. Kays and E. Y. Leung [Ref. 21] 
meported solutions for turbulent flow heat transfer in a con- 
vic Circular tube annulus with a fully developed velocity 
profile and constant heat rate per unit length, for a fluid 
Brandt! number 0.7. Results for large aspect ratio channels 
Meer rectangular cross-sections for laminar flows were reported 
ByeM. Durao [Ref. 22] and J. Ballard [Ref. 23], while R. Holihan, 
Jr. reported results for laminar and transition flows [Ref. 
ШИ Р. Ғ. Brinich апа В. М. Graham [Ref. 25], reported 
results for turbulent flows in a rectangular curved channel 
With an aspect ratio of 6, for the inner wall heated, the 


outer wall heated, and both walls heated. 
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II EN MENO BS TAB STUDY 


The purpose of this investigation was to examine the 
effect of streamwise curvature on the heat transfer rate in 
Tea rectangular duct of large aspect ratio. The flow 
velocities examined were in the turbulent regime. The results 
of this study were compared to the heat transfer rates for 
mE Same range Of velocities, in a straight duct of identical 
aspect ratio. 

Taylor-Gortler vortices are considered to be a laminar 
flow phenomenon. Enhanced heat transfer in curved channels 
for transition and turbulent flows has been observed. Even 
ПИ it has not been proven that the vortices continue to 
propagate at the higher flow velocities, this enhancement in 
the transfer of heat has been attributed to the secondary 
flow velocity components of the Taylor-Gortler vortices. It 
is believed that these secondary components transport the 
heated fluid from the outer concave wall, inward toward the 
opposite wall of the channel, displacing the cooler fluid 
particles and causing them to move toward the heated concave 
wall. It was expected that similar results would be observed 
in this study. 

This investigation was conducted using a single channel 
With a rectangular cross-section and constant aspect ratio. 


The channel incorporated both a straight test section and a 
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curved test section. The results obtained in each test section 
approximately the same flow rates were compared in an effort 
to determine the effects of the Taylor-Gortler vortices on the 
transfer of heat. Also, the results of both the straight and 
curved sections were compared to the results of Ballard [Ref. 
Папа Holihan [Ref. 24]. 

Mees straight section results of this study were compared 
 ГПе results of Shibani and Ozisik [Ref. 19] for turbulent 
flow between parallel plates. In addition, the results of 
Kays and Leung [Ref. 21] and the Dittus-Boelter equation 
[Ref. 26] were compared for Reynolds numbers in the range of 
this study. Kays and Leung's data was for turbulent flow in 
annular passages while the Dittus-Boelter equation was for 
ШОШ Т епс flows in circular tubes. 

The curved section results of this study were compared 
Meee the results of Brinich and Graham [Ref. 25] for turbulent 


flow in a rectangular curved channel. 
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I pe APERIMENTAL WORK 


DESCRIPTION OF THE APPARATUS 

To meet the objectives of this investigation, a channel 
of rectangular cross-section was used. Details of its con- 
EU ton are described in references 22 and 23. Plexiglas, 
0.635 centimeters thick, formed the inner and outer walls of 
the channel. These walls were separated by 0.635 centimeter 
Spacers that also served as the sides of the channel. 

The channel, shown in Figure 5, had a straight section, 
ҰЙ! centimeters in length, followed by a curved section of 
180 degrees of arc, followed by another straight section, 91.4 
E meters in length. The inner curved wall of the channel 
had a radius of curvature of 30.5 centimeters. The rectan- 
mar channel was .635 centimeters high and 25.4 centimeters 
wide resulting in an aspect ration of 40. The hydraulic dia- 
meter was 1.229 centimeters. The cross-sectional area of the 
channel was 16.13 square centimeters. The wetted perimeter 
was 52.07 centimeters. Figure 4 shows a cross-sectional view 
Of both the straight and curved sections of the channel. 

An entrance bell constructed of plexiglas was connected 
to the straight section of the channel. It was designed and 
manufactured according to ASME nozzle standards, with an 
elliptical curved base on a major axis equal to ten inches 


Bud a minor axis of one inch. Cheese-cloth was attached to 
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EN Nutrance nozzle to prevent foreign matter from entering 
the flow channel. An aluminum exhaust nozzle was attached 
ШЕРТПЕ exit of the channel and directed the flow from the 
Siammel into two inch diameter pvc piping. The two inch 
Mvempiping contained a concentric orifice with a diameter 
E065 inches and was constructed according to the ASME 
IE cst Code. Pressure taps on either side of the orifice, 
Band 1/2 pipe diameters from the orifice, and connected 
to manometers were used to calculate the mass flow rate, of 
working fluid, which was air at room temperature. 

The flow of air was drawn through the channel by an 
Сау атіуеп Spencer Turbo Compressor, 50 horsepower 
ПО тртп, гатеа аї 550 cubic feet per minute at 70 degrees 
fahrenheit and one atmosphere. The channel contained two 
test sections from which experimental heat transfer data 
could be obtained. The straight test section was located at 
a distance downstream of the entrance bell that would ensure 
Eusrodynamically fully developed flow. The straight test 
EH was 29.2 centimeters in length with a heated test 
Section area of 741.7 square centimeters. The curved test 
section was located in the lower half of the curved portion 
Or the channel. It was 28.3 centimeters in length, subtend- 
ing an arc of 53.1 degrees, and having a heated test section 
area of 718.4 square centimeters. 

In each of the heated test sections, Temsheet, a carbon 


impregnated porous paper with a uniform electrical resistivity, 
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was glued to the outer wall of each section.  Joulean heating 
sed tO heat the flow of air through the channel. Since 
the electrical resistance of the Temsheet is not constant, 
but varies slowly with temperature, a precision resistor with 
an electrical resistance of 2.0262 ohms was connected in series 
with the Temsheet to allow the calculation of the instantaneous 
power being supplied. Detailed schematics of each test section 
are shown in Figures 5 and 6. 

The variables that were measured and used in this invest- 
ligation were: 

(1) the temperature of the air entering the channel (Ti,) 


) 


E the temperature Of the air leaving the channel (И 

memene temperature Of the heated outer wall for each 
test section Ge) 

(4) the temperature of the unheated inner wall for each 
test section (Te 

(5) the temperature between the three layers of insulation 


) 


(6) the voltage across the precision resistor С 


at each test section GE rs 


(7) the voltage across the Temsheet or heater (Vj) 


(8) the temperature of the air at the orifice (Torig} 


EI) the pressure upstream of the orifice У 


(10) the difference in the upstram and the downstream 
pressure across the orifice (AP) 


) 


(11) the atmospheric pressure (Рат 
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Figure 6. Curved test section, detailed schematic. 
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Б еп о: (пе temperature measurements was made with Copper- 
Constantan thermocouples. Those thermocouples that could be 
removed easily from the channel were calibrated using a 
ROSEMONT Commutating Bridge model 920A, and a ROSEMONT Constant 
Temperature Bath model 913A. A second order polynomial was 
вес tO each set of data. The resulting coefficients dif- 
fered by less than 1.0 percent, so an average of each coeffi- 
Gent was used in converting emf to degrees centigrade in the 
reduction of the data obtained. 

A total of sixty-one Copper-Constantan, 30 gauge thermo- 
Couples were located throughout the test apparatus to record 
the desired temperatures. Four thermocouples were connected 
ШШ рата11е1 to read an average temperature at the entrance 
Симе Е to the channel. Five thermocouples were connected 
in parallel and inserted between each of the layers of insula- 
tion at each of the test sections. In addition, individual 
thermocouples were positioned in each of the test sections, 
eight in direct contact with the Temsheet and eight in the 
Plexiglas of the unheated inner wall. Small diameter holes 
in the plexiglas were drilled to allow the installation of 
the thermocouples. One thermocouple was inserted in the 
two inch pvc pipe downstream of the orifice in accordance 
with the ASME Power Test Code to record the temperature of 
the air flowing through the orifice. The beads of the thermo- 
couples in contact with the Temsheet were electrically insu- 
lated from the Temsheet with ENMAR Heat Resisting Glyceryl 


палате. 





The thermal insulation consisted of 1/4 inch layers of 
Armstrong Armaflex 22 Sheet Insulation, a flexible foamed 
Mastic material, on the heated side of each test section, 
and 1/2 inch Armaflex 22 insulation elsewhere on the test 
apparatus. Three layers of the 1/4 inch insulation were 
used and covered an area slightly larger than each of the 
heated areas. By positioning thermocouples between these 
insulation layers, the heat lost to the environment could 
memeomputed at each test section. The heat loss through 
the 1/2 inch insulation surrounding the rest of the channel 
Че to the heated air flowing through the channel was small 
In comparison. See Figure 7 for a detailed sketch of the 
thermocouple placement and attachment. The insulation was 
held in place by the use of Adhesive Heat Resistant ven- 
tilation Duct Tape. The entire channel and all connections 
between the channel, the pvc piping, and the orifice were 
sealed with General Electrical Silicone Rubber Sealant Caulk 
to ensure there was no leakage of air into the channel or the 
Eu uns which could affect the temperatures recorded as well 
as the mass flow rate. Two aluminum electrodes, 0.318 centi- 
meters thick, were inserted in between the Temsheet and the 
plexiglas in each section as shown in Figure 7. 

The data acquisition system used for this investigation 
Was a Hewlett Packard 3054A Automatic Data Acquisition/Control 
System consisting of a 3456A Digital Voltmeter and a 3497A 


Data Acquisition/Control Unit. Also used in conjunction with 









lemsheet 


\— Copper - Constantan Thermocouples 






Electrically Insulating Glue 


Insulation 





Plexiglas 


1n test sections. 


Thermocouple placement 


Figure 7. 
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BuSOdata acquisition system, were a Hewlett Packard 9326 
Computer terminal and 2671G Printer. Pressure measurements 
were taken from Meriam vertical manometers, one whose fluid 
was water with a 0 to 60 inch range, and the other whose 
fluid was mercury, calibrated to read inches of water with 

a range of 0 to 415 inches. A photograph of the channel and 


meeociated test equipment are shown in Figures 8, 9, and 10. 


EXPERIMENTAL PROCEDURES 

Experiments were conducted by varying the mass flow rate 
of air corresponding to Reynolds numbers from approximately 
12000 to 23000, where the Reynolds numbers are based on hy- 
draulic diameter. For each mass flow rate, data was taken 
automatically by the data acquisition system and the data 
was immediately reduced and printed for examination and com- 
Beeson. The experimental procedures followed were the same 
for the straight and curved test sections. The results ob- 
tained from the straight test section served as a baseline 
for the comparison of the curved section results.  Prelimi- 
nary runs were performed to determine the time required for 
EE west rig to come to steady state. For this experiment 
three hours were allowed for steady state to be reached. 
After three hours, data was taken at ten minute intervals 
to observe changes in parameters to ensure that a steady 
state condition had been achieved. The criteria for steady 


State was based on three variables: 





(1) the mass flow rate 

(2) the heated boundary wall temperature 

E the heater voltage 
When these variables varied by less than two percent over a 
ten minute interval, it was considered that a steady state 
condition had been achieved. 

Mas also necessary to determine the approximate power 
setting required to bring the heated wall of each test section 
to approximately 50 degrees Centigrade. It had been deter- 
mined by Holihan [Ref. 24] that this temperature was sufficient 
to ensure a 20 degree Centigrade difference in the heated wall 
ШК пе unheated wall. 

SS each run, the atmospheric pressure, the pressure dif- 
ence across the orifice, and the pressure upstream of the 
 ШГІсе were entered into the computer for the mass flow rate 
E ulations. А11 other data was acquired by the data acquisi- 
meme system directly from the experimental apparatus. 

To compute the instantaneous power supplied the following 


lationship was used: 


ШИЕ Әтесісіоп resistor voltage and the heater voltage were 
mead by the data acquisition system and the resistance of the 
precision resistor was a known constant. A diagram of this 


emeut is shown in Figure 11. 
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ig SS ERE SENTATION OF DATA 


BE ANALYSIS 

The constant heat flux surface for the heated wall of each 
test section was approximated by the uniform electrical resis- 
tivity of the Temsheet. The insulated unheated wall was con- 
Sidered adiabatic since the heat losses through that wall were 
negligible. The channel's large aspect ratio of 40 allowed 
KE modeling of the channel as parallel plates. The straight 
(топ 0: the channel upstream of the straight test section 
Was of sufficient length to ensure that the flow was hydro- 
 Шшіса!т!ту developed for the flow velocities of this study. 
The straight portion of the channel downstream of the curved 
E UNsectron was also of sufficient length to ensure that the 
flow exiting the channel was thoroughly mixed and that the 


Semperature, T was an average bulk temperature. 


entis 

Based on these assumptions, the experimental configura- 
tion was modeled as forced convection between parallel plates, 
for hydrodynamically developed and thermally developing flows, 
subjected to a constant heat flux. The boundary conditions 
are one wall at a constant heat flux and the other wall 
insulated. 

To analyze this problem, several quantities were defined 
as follows: 


The heat convected to the air was calculated using the 


expression: 
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ee Ср е 


where "С разг. а С пенале те heat of the air at constant 
pressure and 'm' was the mass flow rate of the air. 
The average heat transfer coefficient between the heated 


wall and the flow of air in the channel was defined by the 


eat LON: 
= UE 
И EE 
Арі, AT 
where 'Q,i,' is defined above, 'App' was the area of the 


Temsheet in the test section, and 'AT' was the difference 
between the average heated wall temperature OE o) and the 
average bulk temperature (Ty 14) - The average bulk tempera- 
ture was defined as the arithmetic mean of the entrance and 
J 


The average Nusselt number was then calculated as follows: 


Ew temperatures a ce ios 


h D 
а= Һа 





K_: 
air 


In this equation "ра! is the hydraulic diameter and Е 


Ш the thermal conductivity of air. 
The Reynolds number was calculated for each test run as 


follows: 
м Бра 
А 


Кер 4 = 


c Hair 


40 


where again 'm' and BE are the mass flow rate and hydraulic 
diameter of the channel, pes 1s the cross-section of the 
channel, and И is the dynamic viscosity of the air. 


For the curved section runs a Dean number was defined as: 


Dad 


De = Керч UI 


mere Re' is the Reynolds number based on hydraulic diameter, 
E. Nthe hydraulic diameter, the 'Ri' is the radius of 
curvature of the unheated inner convex wall. 

Ne ten ot the control volume and a set of sample calcu- 
ШЕ ШОП5 for one test run of the curved section are given in 


appendix B. 


Pee RESULTS 

The data obtained from each experimental run was evaluated 
ino the expressions described in Part A above. The major 
parameters resulting from this evaluation are shown in Tables 
ЖЕН ІІ. Table I contains the straight section results and 
EDUC II contains the results of the curved section. A plot 
of the average Nusselt number verses Reynolds number is given 
E scure 12 for the comparison of the results. Error bands 
have been indicated as a result of an uncertainty analysis. 
A sample calculation for the uncertainty analysis is given in 
Appendix A. 

rl Indicate al increase in the rate of heat 


E Ster with increasing Reynolds number for both the 
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ABLE. | 


SUMMARY OF STRAIGHT TEST SECTION RESULTS 


бета Rair 
(W) 
12940 114.43 
14520 I5. 55 
17020 120.00 
20480 ВИ. 57 
22980 qne sl 


Reynolds number and 


diameter. 


11 


(W/m C) 


05.55 


67 


405 


70; 


84. 


2, 


Nusselt 


70 
30 
04 


number 


A Nu 
(C) 
24.27 pM 
24202 Sues 
21.09 36.1 
18.76 BOE 
76 45.4 


are based on hydraulic 





Кеја 
13490 
16950 
#1140 
25520 


ОДА А OFT CURED TEST SECTION RESULTS 


ре 


ЧОТ 


(W) 
С 
50 
114.94 


113.43 


TABLE II 


h 
(W/m C) 
71225 
91.26 
705.246 


111247 


AT 
(C) 
20.79 
17,98 
15, 08 


14.17 


307 
42. 
49. 


Зи 


Reynolds number, Dean number, and Nusselt number are 


based on hydraulic diameter. 
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Nn ssHt' and curved test sections. In addition, the heat 


Bpansfer rate was higher for each Reynolds number investigated 


Еле curved section. 


In the investigation by Holihan [Ref. 24] a straight line 


Meese squares correlation of his data resulted in: 


for 


för 


Nu = 0.983 Re, °°? 300 < Re, < 2000 
laminar flows, and 

Е 036 = 
Nu = 0.110 Ке 4 2000 < Ке « 3000 


transition flows in the straight test section. The curved 


section data yielded: 


BOT 


for 


BOT 


for 


ШІ 0. 065 B 300 < Re, < 2000 
laminar flow, and 

E 0.58 _ 

Nu = 0.117 Re, 2000 < Re, < 3000 


Imansition flows. 


ime data for the present study resulted in 


Nu = 0.063 Re JE 10° < Ке, а < 2.5 x 10° 


hd 
turbulent flow in the straight test section, and 


т- _ ШЕСІ 4 
Мп - 0.040 Ке, 4 10 < Ке, 4 с: < 10 


4 


star flow in the curved section. 
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In the turbulent regime of this study, the rate of heat 
Master was about twenty percent higher in the curved section 
Dueuscompared with the straight section. This compares favor- 
ably with earlier studies by Kreith [Ref. 16], where he re- 
Mertead an increase in the rate of heat transfer along а 
 псауе wall of from twenty-five to sixty percent for 
Reynolds numbers, based on hydraulic diameter, between 


| and 109. BINE асро an increase of fifteen per- 


10 
E or laminar flow and thirty percent for transition 
flow [Ref. 24]. Ballard reported an increase of eleven 


percent in the heat transfer rate for laminar flows [Ref. 23]. 
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V. DISCUSSION AND CONCLUSIONS 


Holihan [Ref. 24], in his study, determined that there 
was a negligible difference in the fluid bulk temperature 
and the unheated wall temperature. He additionally demon- 
Strated that the radiated heat transfer from the unheated 
Wall surface was minimal. These results were verified for 
МИР ртезепіт study and provide the basis for the assumption 
Maat the heat transfer to the air flowing through the channel 
was solely by convection from the heated Temsheet. 

As was mentioned earlier, the high aspect ratio of the 
@aammel provided the basis for the assumption that the ex- 
mental apparatus, as configured, could be modeled as 
infinite parallel plates. Holihan's data [Ref. 24] also 
tended to verify this assumption. His experimental data, 
йе laminar flow region, approached the theoretical limit 
for average Nusselt number for parallel plates with one wall 
Meated at a constant heat flux and the opposite wall 
Matabatic [Ref. 27]. 

Based on the assumptions mentioned above, comparisons were 
made with other analytical and experimental results that were 
of the same or similar problem. Again the problem being, flow 
between infinite parallel plates with one wall at a constant 


heat flux and the opposite wall adiabatic. 
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comparison of the experimental results of this study 
and the experimental results of Ballard [Ref. 23] and Holihan 
[Ref. 24] are shown in Figure 15, for the straight section, 
Нисте 14 for the curved section. Figure 15 shows a com- 
Ata tion ot all the data on one plot. Ballard's data was 
ШОСТУ іп the laminar region, while Holihan's data covered 
the laminar as well as transition regions. The channel aspect 
ratio was 40 for each of these studies and the experimental 
cedures were similar in each case. Ballard and Holihan 
ррпратесй their data to the analytical studies by McCuen, 
[Ref. 28], for heat transfer between infinite parallel plates 
with constant wall temperatures and heat flux. Their results 
plotted above the analytical solution and the difference was 
attributed to the difference in the geometries used in the 
вене, and to the limitations inherent in any experimental 
E Ihe side wall effects and the inability to totally 
mount for all the heat transfer processes and/or losses 
were also factors that may have caused the difference. 

A comparison of analytical and experimental results with 
Base of the present study for the straight test section is 
shown graphically in Figure 16. The results of Shibani and 
Ozisik [Ref. 19] is shown for turbulent flow between parallel 


plates. The correlation used was 


ee 024 ве 82 
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IU! < Pr « 10 and 10° < Re < 109. The Dittus-Boelter 


equation [Ref. 26] for heat transfer in a straight tube for 
constant wall temperature is plotted from the following 
@orrelation: 


Nu = .02 Re^? 


for a Prandtl number equal to 0.71. Additionally, Kays and 
Leung [Ref. 21] obtained data for heat transfer in annular 
passages, with one wall heated at a constant heat rate. 

The data shown is for r* = 1.0 which equates to parallel 
puates. 

Nb curve section data is plotted in Figure 17 and is 
compared with the experimental results of Brinich and Graham 
[Ref. 25] for flow and heat transfer in a curved channel. The 
accuracy of the data points used are subject to errors, in 
that the actual values were not given in their study and had 
вое taken from a plot of Stanton number versus Reynolds 
number as given in the reference. 

The data fromthe present study fell below the correlations 
that were used as comparisons. These differences can also be 
Eu buted to the differences in geometries and conditions of 
Mating used. Experimental error and side wall effects may 
Meer aiso contributed to the overall difference in results. 

Ше Пата from the present study and that from the Brinich 
and Graham study seem to compare favorably and exhibit a simi- 


IT trend. Irrespective of the differences noted above, it 
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Ts that there is a definite enhancement of the heat 
Branster rate in the curved test section for laminar, tran- 
Sition, and turbulent flows. Streamwise curvature, result- 
Bun secondary flows, seems to cause the enhancement. 
Whether the Taylor-Gortler vortices continue into the 
transition and turbulent flow regimes remains to be seen. 
Continued study of the effect of streamwise curvature and 
meemaevelopment and propagation of the Taylor-Gortler vor- 


Не throughout all flow regimes is required. 


55 





VI. RECOMMENDATIONS 


impetter understand the effects of streamwise curvature 
meee enhancement of heat transfer, a significant amount of 
ШИ 50111 remains. With the introduction of today's data 
Meqmisition systems, in conjunction with portable computers, 
Instantaneous recording of data and its reduction, for imme - 
diate analysis, saves an immeasurable amount of time and 
meamces the chance of human error. The present experimental 
apparatus has been used over a range of Reynolds numbers 
throughout the laminar and transition regimes, into the tur- 
bulent regime. Results have compared well with previous data 
E Correlations. A new flow channel of similar design should 
Memeonstructed prior to further testing, due to the age and 
Вот стоп of the present equipment. A new set of calibrated 
thermocouples and possibly the introduction of pressure trans- 
ducers, instead of manometers, that could be connected to the 
data acquisition system, would provide more consistent, and 
as a result, more accurate data for flow rate measurements. 
Also, since the present apparatus has been modified several 
Mees in the past, resulting in leaks that are difficult to 
locate and stop, a new flow channel would help to eliminate 
Some Of the inaccuracies to temperature and flow measurements 
introduced by undiscovered leaks. It is also highly recom- 


mended that a method of safeguarding each piece of equipment, 
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вас in the event of a power failure, when power is re- 
(Әй, each item would require the operator to manually 


restart it. 
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КОРГЕШ ОТКА: EXPERIMENTAL UNCERTAINTY 


Tm uncertainties for the major variable in the experiments, 
meme, Calculated in accordance with the method described by S. 
Mime and F. McClintoch [Ref. 29]. The estimates of the un- 
certainty in the measured quantities were made conservatively. 
As a result, there is considerable confidence in the uncer- 
tainties as calculated. 


The following equations were used to calculate the 








uncertainties: 
17/2 
7 2 à 
И сл - а [ука (5) + (40 + (50 ) 
2 Ае 
+ (1/2 0,; AP T ашу, ES 1 
d 
2 2 IT 
_ dm air 1 out “in 
ПАО SQ, (O + Mv LÁ 
1/2 
dQ... 2 dA,, 2 
(3) dh» E (qu + a ат) 
ЗЛ PL 
1/72 
2 dD dk =. 2 
(а) амы = Ка (098) СТВО + (1880) 
h Did au 
172 
D d 2 
hd acer G 
БЫС -5 (23M + + E Se 
ES ) E ( E. Á 
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The values obtained for the uncertainties for the values 


in the curved section test run at a Reynolds number of 


8500. аге: 
Quantity Uncertainty 
A ‚0005 
Арі, SOO AG 
A. .0098 
Cp .0054 
Did 0087] 
h .0412 
K .0050 
Ir 0012 
Nu .0421 
m 0225 
Е... „0582 
Rey gq .0260 
Thik 2:019 97 
n .0046 
Bu .0049 
IE -0065 
в т... „0150 
ү -0051 
AP 20227 
AT 0057 
ШО ШЕ 
Ee. .0005 
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ihe Major source of uncertainty in the average Nusselt number 
is the uncertainty in the pressure readings from the mano- 
meter which affects the mass flow rate calculation. In 

order to obtain a more accurate pressure measurement, it 
eocen recommended that pressure transducers or inclined 
manometers be used instead of vertical manometers in future 


studies. 
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NEE DIN EB SAMPLE CALCULATIONS 


Figure 18 below, shows the major heat transfer components 
mS eel of the test sections. The sample calculations that 
ПОШТОМ, demonstrate the methods used by the computer to calcu- 
late these components as well as the Reynold's number, average 
heat transfer coefficient, and average Nusselt number for each 
ET of data. The sample calculations are for the curved sec- 
tion, but those for the straight section are similar. Reynolds 
number, Nusselt number, and Dean number are based on hydraulic 


diameter. 





o" Qu 


РЕ Бпегоу balance in straight section. 
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SAMPLE CALCULATION DATA 
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Patm 
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E 
3 
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He. 


insl 


ins2 


pj и О 


Grif 
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air 
Chair 
air 
Kins 


ans 


gc 
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35% 


4. 


2 
50 


о 
1-55 
7 зм 
28. 
44. 
гою 


Son 


о 


28% 


26 


92 


To s SV 


4075 V 


20202 10 


ӘРІП He 


sein. H,O 
NS. H,O 
41 
08 
50 
99 


66 


(C Же COSME) Mu оо 


52 


аак 


1006 J/Kg К 


18 


SE 


d 


1 


roo ox 107 


ECCE Г Кол 


6 Kg/ms 


2 


6495 x 10 ^ W/mK 


ONS Mm 
70 


. 40 


4 


23925 


402 


5115. 


286.8 Nm/Kg K 
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Mo vi = 1.0 
4 2 
Apr, = ‚07188 т 
А = .0016 п 
C 
n = .002027 m? 
pipe 
D. = .00635 m 
D. = ,0508 m 
pipe 
Пе = ОЕ п 
ot = xm 


TEMPERATURE CALCULATIONS 


1. Bulk Temperature (Ty 14) 


2. Mean Temperature Difference (AT) 
AT = Db 2014500207752 17.75 C 


) 


5. Temperature Difference in the Insulation (AT;,. 


AT. ни SE 


= = x 2 - 
ins ШЕТ 38200 252 ОО ТАС 


IRS 2 


POWER CALCULATIONS 


1. Power Supplied (Q5) 
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ШЕН 5509525904.4075) _ 


120.14 W 
P RPR (2.0262) 


cat Lost Through the Outer Plate (Q,,) 


| АШЫ е 6.14 
= > Бе о О пола полом 
ШТА ARK Ay (0.00635)/(3.6495x10- 4) (.07188) 
- 2.54 W 


3. Heat Radiated (059) 


a. Radiation Resistance (Вр) 








T. 1 “был | 1 ] p 
|o ل‎ 
R ApLEwO ке -wi "Арі wi 
T |: 1 1 ] ] 
= nu [ + = 1] = туртту ото 1 = 40.74 
PL “wo Ewi у ° | -2 


b. Heat Radiated coe) 


(5.669x10"%) [(44.50+273)* (28 
ee ——— 


IZ 


D. MASS FLOW RATE CALCULATIONS 


1. Pressure Conversions 


а. P,.4,730.45 in Hg x 3374.1 = 102741.35 N/m* 
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E. 


p 


_ 
БОО ОАР = 22.5 in H50 x 248.64 = 5594.4 N/m 


g 
І 


с E ШЕН 218 54) 4 102741.55 - 105980.10 
2 
N/m 
Density wot Air De) 
_ P atm E EI AS) NL E 3 
air PLANT TE КЕЛЕР . AS m 


Grit 


Expansion Factor (Y) 


ү = EINE 1450930. 78-128”) ] aD 
1 


| 


E TI 1T М0 7(.5325)?+12(.5325)-°] 


5594.4 
X (1.402) (105986.10) 


0758 


Area of Orifice (A) 


E. "orig ERU — sonas 2 
4 4 
Mass Flow Rate (m) 
ii = YKA YZgep,,, BP = (.9738)K(.0005747) 
ЕЙ ЕТІ RSS Sr 
й = .06455 К 
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il crating: 
Assume a Reynolds number pipe - 55000 
Obtain a value for K, the flow coefficient, from 


reference 30. K = .6332 
Solve for m. m = .0409 Kg/s 


aD. 
Solve for new Re.. = П іре = (.0409) (.0508) 
Ds 1.002027) (18.63х10 9 


pipe air 


R = 54988 


“pipe 
Check convergence and repeat process if necessary. 
(Convergence 1f difference less than .001) 


m = .0409 Kg/s 


REYNOLDS NUMBER CALCULATIONS 


h D. 
1. Regg, = RES - UgUZUTH IS.  xIU-6 7 54988 
БЕР pipe ЗЕЕ Й 
m D 
E. с _ 0.0409) (.00635)  _ g713 
а Ас Hair . 1 1o OSX 
m D 
hd 
3. Rega * FT. 
С arr 
4 x< A 
а. бал کس‎ = A + 201229 n 
wet 
M s Dhd _ (.0409)(.01229) Т 
һа DNE. (.0016)(18.63x10-6 


ШЕУ CONVECTED TO AIR CALCULATION 


1. Фају = АС (Тоџе Tey) 7 (.0409) (1006) (28.08-25.41) -109.86W 
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СЕБЕ HEAT TRANSFER COEFFICIENT CALCULATION 


Q В 
A= air . 109.86 . $6.11 W/m^C 


Арт АТ ШИВ (17.75) 


AVERAGE NUSSELT NUMBER CALCULATION 





h D 
NT = е ШЕ p - 59.99 
ee 205927 10 


DEAN NUMBER CALCULATION 


Dna 


E р 01229 _ __ 
8 = 16863 Feet 5445595 


Пе - Ве 





һа 
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APPENDIX C: CORRELATIONS 


The correlations obtained for the present study were 
obtained using the method of least squares for a first degree 
poesmnomial, as outlined by C. Е. Gerald [Ref. 29]. In this 
method the values of Reynolds number (x) and Nusselt number 
(Y) were first converted to their natural logarithmic value. 
2 


Next the summations of NES 2 Jac and XY. were calculated. 


The values of these quantities were then placed in a matrix: 


a ЕЭ 
N =x. О d 
Me 21, 1=[ | 
хі їх; jr DX тү 
and the resulting simultaneous equations 
ав + a, Xs С 
2 


а. 2Х1 + a] 2X; - LX ¿Y 


were solved for a. and ај. 


The equation of the line then became: 


Ln Nu = a + a, Ln Re 
Or 


№ = ео ВЕ | 


Mme standard error for the present study correlation was 


Calculated by the following equation: 


М - п - 1 


С = 
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Where: Yi is the actual value of the average Nusselt number 
obtained experimentally. 
y; is the value of the average Nusselt number calculated 
from the correlation equation. 
Ес the number Of data points used in the correlation. 
nis the degree of the polynomial used. 
oS the present study the following correlation was obtained 


EE straight section: 


— 0.65 
Nu = ‚0635 Кеја 


ENS standard error of: 


5- 0.1988 


The correlation obtained for the curved section was: 


== _ TZ 
Nu що 040 Кера 


and its standard error was determined to be: 


EL 0.1050 
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